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Abstract 

The contact conductance between graphene and two quantum wires which serve as the leads to 
connect graphene and electron reservoirs is theoretically studied. Our investigation indicates that 
the contact conductance depends sensitively on the graphene-lead coupling configuration. When 
each quantum wire couples solely to one carbon atom, the contact conductance vanishes at the 
Dirac point if the two carbon atoms coupling to the two leads belong to the same sublattice of 
graphene. We find that such a feature arises from the chirality of the Dirac electron in graphene. 
Such a chirality associated with conductance zero disappears when a quantum wire couples to 
multiple carbon atoms. The general result irrelevant to the coupling configuration is that the 
contact conductance decays rapidly with the increase of the distance between the two leads. In 
addition, in the weak graphene-lead coupling limit, when the distance between the two leads is 
much larger than the size of the graphene-lead contact areas and the incident electron energy is 
close to the Dirac point, the contact conductance is proportional to the square of the product of the 
two graphene-lead contact areas, and inversely proportional to the square of the distance between 
the two leads. 
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1. Introduction 



The fabrication of graphene, a single atomic layer of graphite, is the first experimental 
realization of realistic two-dimensional crystal In the effective mass approximation, va- 
lence electrons in such a carbon material obey the massless relativistic Dirac equation, rather 
than the Schrodinger equation. Consequently, graphene presents many unusual electronic 



transport properties, such as the half-integer quantum Hall effect 
conductivity minimum at vanishing carrier concentration js, Q, [J oi 



subtle weak localization 1^, [isl . 



3, HQ, 
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6|, the nonzero 



m QJ, the 



carrier through an arbitrarily high barrier [20|, 



and the reflectionless transmission of the 
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The intriguing transport properties mentioned above are closely associated with the scat- 
tering of a Dirac electron by impurities, defects or gated barriers, etc. However, apart from 
these scattering effects, the contact between graphene and the metallic electrodes also in- 
fluences the electronic transport spectrum to some extent. In other words, the measured 
conductance(conductivity) spectrum often includes a contact conductance, accompanying 
the scattering conductance. From the experimental viewpoint, such a contact conductance 
makes sense when graphene couples to the probes of the scanning tunneling microscope 



(STM) 
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25| . or ultrathin gold or tungsten wires; in particular, the multi-probe STM 



that was developed quite recently [26 
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30|, which can be used to explore the 



surface structure of graphene. In these actual structures, the leads that connect graphene 
and electron reservoirs are the quantum wires which have only a few transverse modes to 
carry the current. This is in contrast with the infinitely many transverse modes in graphene. 
Thus, the electron will be reflected with certain probability when it enters the leads from 
graphene. Besides, such quantum wires couple locally to a finite number of the carbon atoms 
of graphene. As a result, the electronic tunneling between the leads and graphene depends 
sensitively on the coupling configuration at the contact. To sum up these facts, it can be 
anticipated that the coupling between graphene and a quantum wire will give a nontrivial 
contact conductance. 

So far, relatively fewer works involve the contact conductance between graphene and the 



leads of different kinds 
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321], in comparison with the scattering conductance. In fact, many 



aspects concerning this issue deserve further explorations. In the present work, we carry out 
a systematic investigation on the contact conductance between graphene and quantum wires 
which serve as the leads to conduct the current. There are many factors that influence the 
contact conductance, for example, the bandwidth and the band position of the quantum 
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wire(lead) relative to the Dirac point of graphene, the couphng configuration between the 
leads and graphene, and the distance and the relative orientation between the two leads. To 
obtain a quantitative dependence of the contact conductance on these factors it is desirable 
to analyse the observed conductance spectrum in relevant experimental setups. 

The rest of this paper is organized as follows: In Sec. 2, the theoretical framework to 
formulate the contact conductance is elucidated. In Sec. 3, the numerical calculations on the 
contact conductance are performed. Based on these numerical results, the dependence of 
the contact conductance on the relevant parameters is discussed. Finally, in Sec.4, we make 
some concluding remarks. 

2. Model and theory 

As illustrated in Fig.l, the electronic transport structure we consider consists of a 
graphene monolayer coupling to two quantum wires, which serve as two leads to conduct 
the current. We adopt a semi-infinite linear lattice to describe the band structure of a single 
transverse mode in each lead. Such a single mode lead couples to the finite number of the 
carbon atoms in graphene. Note: Some structure parameters are explained in Fig.l. 

In such a model the tight-binding Hamiltonian of an electron is composed of three parts, 

H^Hc + Hg + Ht. (1) 

The first term Hq is the Hamiltonian of the two semi-infinite linear lattices which model the 
two leads. It is given by 

Hc^Yl ^^^h + Yl (^oc]+iCj + H.c.) + ^ocjcj- + + (2) 

j<-i i<-i j>+i j>i 

where (cj) is the electron creation (annihilation) operator at the individual lattice points. 
For the ideal leads, both the on-site energy and the hopping energy are viewed as uniform 
parameters, denoted by and vq, respectively. The second term Hq is the tight-binding 
Hamiltonian of the tt band electron in graphene. It takes a form as 

HG^tJ2idlAd,B + }i.c.), (3) 

where dj^^ {diu) with z/ = A or B is the electron crcation(annihilation) operator associated 
with the local atomic orbits in graphene. The notation means that the summation is 
restricted between the pairs of the nearest neighbor carbon atoms, t is the corresponding 
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hopping energy. We choose the Dirac point as the energy reference point, hence the on-site 
energy of all lattice points in graphene is zero. Moreover, in what follows we use the units 
h = t = a = 1. The last term Hj' in the total Hamiltonian describes the electron tunneling 
between the leads and graphene. It is given by 

diy + H.c.) + ^{vajycldjy + H.c), (4) 

iv jv 

where v^iy and v^jy denote the coupling strength between two leads and the individual 
carbon atoms in graphene. Note that we use L and R to denote the two leads respectively. 
As illustrated in Fig.l, only finite carbon atoms around the tip of each lead couple to the 
lead effectively. 

To study the electronic transport properties, we need to calculate the linear conductance 
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34| . we can write the 



between the two leads. Based on the Landauer-Biittiker theory 
linear conductance formula as 

Q{E) = ^Tr[r^(E)G^(E)r«(E)G'^(£;)], (5) 

where E is the energy of the incident electron from one lead. The matrices of the retarded 
and advanced Green functions satisfy a relationship G^{Ey = G°'{E). The matrix ele- 
ment is given by [G%y,jy^ = J^^G^,y^jy^{t)e'^*dt and G^^jVal^) = -i^it){{diu^{t),d]^J), 
just following the standard definition of a retarded Green function. We will often use an 
alternative notation {{AlB))" to denote the retarded Green function in Fourier space for 
convenience, e.g. [G^]iy-^jy2 = {{diuild^juj)^ ■ In Eq.(l5]) two other matrices are defined as 
[T%u,i'u' = -2vLiuVli,^,lm{{c-i\c^_-i))Q and [T%u,j'u' = -2vRjuV*j^j,^,lm{{ci\c\))Q, where the 
subscript "0" denotes that the corresponding Green function belongs to an individual lead, 
completely isolated from graphene. Besides, as we have done, we will often drop the energy 
dependence of these matrices to keep expressions compact. 

To solve the Green functions in the conductance formula, we need to convert the Hamil- 
tonian Hg into the eigen-representation of the vr band electron of graphene. Doing so, we 
utilize the following transformations 

diu = :^$^e*'=-^-Cfc„ (6) 

and 



CkA = J ^fc-B = ^^0(ks, (7) 
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where N is the total number of unit cells in graphene; tk — Yl^=i ^ exp(ifeT;) with k being the 
electron wave vector measured from the center of the Brillouin zone of graphene. The creation 
(annihilation) operator c\j^{cku) is associated with electron Bloch state in one sublattice of 
graphene. And the notation s = + or — , denotes the conduction or valence bands of 
graphene, respectively. al^{aks) is the creation (annihilation) operator of the eigen-states of 
the valence electron in graphene. Thus, the Hamiltonian of graphene becomes diagonal 

Hg = ^Eksal^aks, (8) 

k,s 

where Sks = s\tk\- Accompanying such a representation transformation, the tunneling Hamil- 
tonian Ht can be rewritten as 

Ht = J2i^Lkscliaks + H.c.) + J^^'^RkscWs + H.c), (9) 

ks ks 

with 

Meanwhile, the conductance expression changes into 

2 

4Im((c_i|cLi))SIm((ci|cI))S. (11) 
The retarded Green function satisfies the following equation of motion 

iaks\al,^,)Y = gkshs,k's' + fi'fevVLfe's'fc't^L\s((c-i|cli))'' + 5'fcvVLfc's'fcV'4^((ci|cLi))^ 

+gk's'VRk>s'gksVlks{{c-i\c\)Y + 9k's'VRk's'gksV^Xci\c\))\ (12) 

where Qks = {E~^ —£ks)~^ and = E + irj with rj being a positive infinitesimal. Here, we can 
see that the Green functions {{cxlc^^))^ i^,^' = =tl) ai'c the relevant quantities to calculate 
the linear conductance, which are exactly solvable. Taking ((c_i|c^))'' as an example, the 
set of the equation-of-motion with regard to this Green function consists of 

(£;+ _ - CLL){{c-i\cl,)r = 1 + Vo{{c.2\cU)r + iLR{{c,\cU)y, (13) 

and 

(£;+-£o)((c-(,-+i)|cLi))^ = ^o((c-(,-+2)|cLi))^ + ^o((^^^^ for J = 1,2,.... (14) 

In the above equations 

W = ^Vt^ksgksV*>ks, II, li' or R. (15) 

ks 
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J2 Yl ^Lks{{aks\ak's')y^Rk's' 

ks k's' 
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By virtue of the property of a tri-diagonal matrix 351], we can solve analytically the above 
set of equations, which yields 

with 

uj = E-e^-l3\ (17) 

and 

(5* = E-eo-i.JH-iE-e,Y ^ ^^^^ 

In similar manner, we can obtain the analytical forms of Green functions ((c_i|cj))'', ((ci|cL^))^ 
and ((ci|cj))'". They are given by 

= ... . w./^^ ^ . . . (19) 
= . ^ . . ^ (20) 

and 

ici\c\)Y = i — , r\ , , • (21) 

- Ul>[uj -^rr) - Ur^^rl 
Substituting these results into Eq. (fT2!) . we can then get the explicit form of the Green func- 
tion {{aks\oil,g,)y' ■ Subsequently, after some derivation, we obtain an analytical expression 
for the contact conductance. That is 

g{E) = ^T(E), (22) 
TiE) = \poiE)iiE)\^ (23) 

with 

m = - — , , , (24) 

and 

Po{E) = -2Im((c_,|ct = ^^-'^-{E-eo)\ ^^5) 

So far we have obtained a conductance expression in terms of the self-energy terms ^^^r. 
We will derive their the analytical forms in the Appendix. To calculate the conductance 
spectrum, the two parameters vq and Eq ought to be taken the appropriate values to guarantee 
the energy band of the leads and the linear region of the vr band of graphene overlaps each 
other. Besides, the incident electron energy should also be restricted in the linear region 
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of graphene band. Now we consider an extreme case that the graphene-lead couphng is far 
much smaller than the bandwidth of the leads, i.e. v <^ vq. From the above conductance 
expression, we can infer that the conductance formula in such a weak coupling limit reduces 
to 

G{E) = ^\p,{E)U\'. (26) 

3. Results and discussion 

After formulating the linear conductance, we are now in a position to perform the numer- 
ical calculation for the conductance spectrum, from which we can investigate the electronic 
transport properties dominated by the contacts between graphene and two quantum wires. 

First of all, we consider the simplest case that only one carbon atom in graphene couples 
to each lead(Hereafter we call the tip of a lead a probe). For the numerical calculation, we 
fix the first probe(L) to couple to the A atom at origin, and shift the second probe around. 
In FiglH we show the calculated contact conductance spectrum {G{E) versus E) for the 
second probe (R) stopping at a specific A atom. In such a case we have f^^j^ = v6io6^A and 
VrJu = v5jm5yA- We refer to such a configuration between graphene and two leads as A- 
A coupling. The most striking feature shown in this figure is that the contact conductance 
vanishes when the incident electron energy is aligned with the Dirac point of graphene. From 
the analytical results given in the Appendix, we can readily obtain the self-energy terms ^^^i 
for the case of the simple A-A coupling 

_ y-^v'E qdq 

27r io E+^-{^qr ^ ' 

V3^2^cos(^) n^Jo{qK}qdq_ 

Owing to the chirality of the vr band electron in graphene, these self-energies vanish at Dirac 
point (-E = 0). Then, this leads to the zero point of the contact conductance as shown in 
Figj2l Therefore, we can attribute the zero contact conductance at the Dirac point to the 
chirality of the Dirac electron. As shown in Figl2](a), we can shift the position of the band 
bottom of the leads relative to the Dirac point by varying the parameter sq. As a result, the 
conductance spectrum becomes asymmetric with respect to the Dirac point. In addition, 
Figj2](c) shows that when the band of the leads is widened by increasing the parameter vq, 
the conductance becomes notably smaller. Such a dependence of the contact conductance 
spectrum on the band structure of the leads can be readily understood by analyzing the 



local density of states of the electron at the tip of a lead, which appears in the conductance 
formula, see Eqs. (|2^ . (12^ and From Eq. (12B]l we can see that such a local density 

of states decreases as the bandwidth of a lead gets larger. Besides, when the parameter Eq 
deviates from the Dirac point, the local density of states becomes asymmetric relative to the 
Dirac point. The variation of the conductance spectra with the parameters Vq and Eq shown 
in Figj2](a) and (b) just reflects these features of the local density of states. In Figl2](c) and 
(d) we show the contact conductance spectrum for the so-called A-B coupling which means 
that the second probe couples solely to a S atom at a specific position. In such a case the 
band structure of the leads influences the conductance spectrum in the same way as the case 
of A-A coupling, namely, the conductance spectrum becomes asymmetric with the shift of 
Eo and decreases globally with the increase of vq- However, for the case of A-B coupling a 
conductance peak occurs at Dirac point in place of the conductance zero in A-A coupling. 
This is due to that the self-energy ^^j, for A-B coupling takes a different form. It is given by 

_ V37t^^sin(^ + g) J,{qR)q'dq 

By a simple derivation, we can further deduce that for the A-B coupling the self energies at 
the Dirac point are = = and 

U = ^ ^[1 - MQcR)]. (30) 

Then the contact conductance at the Dirac point takes a simple form as 

GiE = 0) = (4^)^- (31) 

With the help of these self-energy terms we can discuss the dependence of the contact 
conductance on the distance between the two probes. In fact, from Eqs. (!29l) - (!3T!) we can 
infer that initially with the increase of the distance between the two probes, the contact 
conductance near the Dirac point will decay rapidly. However, when the distance between 
the two probes becomes sufficiently large, the contact conductance tends to be inversely 
proportional to R^. In FigJ31(a) we plot the contact conductance as a function of R by 
letting the second probe to move away from the origin (the position of the first probe) 
along the y axis. We can see that the contact conductance decreases drastically when R 
increases within several times the lattice constant. Such a rapid decay has little to do with 
the variation of the incident electron energy. 

The explicit expression of the contact conductance at the Dirac point shown by Eq. fl3ip 
indicates a conductance peak occurs when vq = C,lr- This implies that the resonance will 



occur when the distance between the two probes takes an appropriate value, which depends 
on the bandwidths of the leads and their coupling strengths to graphene. To study such 
a resonance in some details, we calculate the conductance as a function of the strength of 
the graphene-lead coupling, as shown in FigJH From FigHl^a) we can see that a stronger 
graphene-lead coupling is required to observe resonance when the distance between the two 
probes gets larger. The calculated results shown in FigHl^b) indicate that the resonant 
conductance peak is notably suppressed when the incident electron energy deviates from the 
Dirac point. 

Now we turn to discuss the orientation dependence of the linear conductance when the 
second probe shifts around the origin where the first probe is located. For simplicity, we only 
consider the weak graphene-lead coupling limit, where the conductance expression takes a 
relatively simple form, given by Eq. fl26!) . From Eqs. fl28l) and (!29l) we can see that the factors 
cos(^^^) for A-A coupling and sin(^^^ + 6) for A-B coupling determine the orientation 
dependence of the linear conductance. Our derivation in the Appendix indicates that such 
orientation factors originate from quantum interference between the K and K' valleys. Con- 
sequently, when the length of -R(the distance between the two probes) is fixed, the maximum 
of the conductance appears when R is along y direction. Furthermore, considering the sym- 
metry of the honeycomb lattice of graphene, when the second probe moves along a circle 
around the first probe, the conductance maximum will appear at the ir^ directions. In 
addition, we can readily infer that for A-A coupling the conductance pattern formed by 
shifting the second probe has reflection symmetry with respect to the x axis. In other words, 
when the second probe is located at the two symmetric A atoms relative to the x axis, the 
conductance gives the same value. 

The above coupling manner that a quantum wire couples to a single carbon atom in 
graphene can only model the extreme situation that the probe aims at a specific carbon atom, 
but the coupling strength between the probe and the adjacent carbon atoms is negligibly 
small. To mimic the actual graphene-lead coupling configuration, we introduce the following 
Gaussian-type graphene-lead coupling 



where r denotes the arbitrary position of probe /i, and Rii, is the lattice vector of a carbon 
atom in the vicinity of the probe, do is a decaying factor to determine the region in which 
the carbon atoms couple effectively to the probe. Obviously, when is far smaller than 
the lattice constant, such a Gaussian- type coupling changes into the single atom coupling 




'0 



(32) 
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configuration discussed above. On the other hand, when do is comparable to or larger than 
the lattice constant, a probe will couple to multiple carbon atoms around it. 

At first, we check whether the contact conductance vanishes at the Dirac point when the 
two probes are positioned at two carbon atoms of the same kind. As shown in Fig|3](b), we can 
see that such a chirality associated feature no longer exists when the effective coupling region 
gets larger. This is readily understood since in the Gaussian-type coupling configuration with 
a decay factor do comparable to the lattice constant of graphene, the conductance is the 
averaging result of A-A and A-B couplings. Besides, the result shown in Figl3](b) indicates 
that the rapid decay of the contact conductance with the increase of the probe interval still 
holds in such a Gaussian-type graphene-lead coupling. In FigJl](c) we show the resonant 
contact conductance by changing the strength of the graphene-lead coupling in Gaussian 
type. We can infer that the resonant conductance peak is influenced by the decay factor do 
intricately. Initially when we increase the decay factor((io = 0.4 to 0.8) there is a notable 
decrease in the resonant peak. But when we increase do further, the conductance peak turns 
to get larger. Such a complicated dependence of the conductance on the effective coupling 
size arises from the quantum interference between the different electron transmission paths. 
When an electron travels from the first probe to the second one, electronic tunnelings between 
graphene and the two probes can be realized via two carbon atoms of either the same kind 
or the distinct kind, which corresponds to the different electron transmission paths. These 
two distinct electron transmission paths result in destructive interference. At do = 0.8, 
besides the central B atom, three nearest neighbor A atoms begin to couple effectively to 
the second probe. As a result, the destructive quantum interference diminishes the resonant 
peak. When do = 1.2, six next-nearest neighbor B atoms enter the effective region to couple 
to the second probe, which compensate for the negative effect of the three nearest neighbor 
A atoms. As a result, there is a notable increase in the resonant peak. 

With such a Gaussian-type coupling manner we can plot the two-dimensional conduc- 
tance pattern formed by shifting the second probe around the first one. In Figj5] we plot 
such a conductance pattern for a relatively small decay factor in the Gaussian-type coupling 
function((io = 0.2). This case is in analogy to the single atom coupling configuration dis- 
cussed above. From this figure we can see that the conductance pattern exhibits the C3 
group symmetry when the second probe rotates around the first probe. And such a rota- 
tion symmetry is independent of the incident electron energy. The strongest conductance 
appears when the second probe is located at the next-nearest neighbor B atoms, rather than 
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the three nearest neighbor B atoms around the origin (Except for the conductance maximum 
at the origin when the incident electron energy is far away from the Dirac point, as shown 
in Figin](b)). According to the discussion on the relation between the resonant conductance 
peak and the probe interval, the strongest conductance can appear at different positions if we 
adjust the coupling strength between the probes and graphene. In Fig|6]we plot the conduc- 
tance pattern corresponding to a relatively large decay factor((io = 1-0). We can see that the 
C3 group symmetry remains in this conductance pattern. And conductance pattern does not 
vary notably with the variation of the incident electron energy. Unlike the case of do = 0.2 
where the conductance maxima always occur when the second probe points at an individual 
carbon atom, some discrete islands with different shapes form in FiglHl labeling the regions 
with the relatively large conductance. As shown in Fig|7l if we move the first probe to the 
center of a hexagon of the graphene lattice, the conductance pattern formed by shifting the 
second probe around has the Cq group symmetry instead of the C3 group symmetry in the 
previous cases shown in FigjS] and FigJHl However, in such a case the conductance is far 
smaller than those in the previous cases. This is because that the first probe does not point 
at any atom carbon. Thus, according to the Gaussian-type coupling function, the coupling 
strength between the carbon atoms in graphene and this lead diminishes notably. 

Finally, based on the above discussion, we can work out an approximate expression for the 
contact conductance, which is convenient for a rapid estimation of the contact conductance 
for some experimental setups if the following conditions are satisfied. First, if the distance 
between the two probes, symbolized by R, is much larger than the size of the effective 
contact area around each probe, thus to calculate the self-energies from the formulae given 
in the Appendix, we can replace the exact distance between two individual carbon atoms 
around two probes by R approximately. Second, when the incident electron energy is limited 
in the vicinity of the Dirac point, we can roughly view E ^ 0. Thus, according to the 
formulae shown in the Appendix, the dominant contribution to the self-energy arises from 
F^g^R) and FgA^R). Consequently, we arrive at an approximate expression about the 
contact conductance in the weak coupling limit, which is given by 

2e^ n']niv'^ 



h Avl 



\Fba\\ (33) 



where denotes the number of the same kind in graphene which couple effectively to probe 
/i. Obviously, such a quantity is proportional to the contact area between graphene and the 
probe. If the probe interval R is much larger than the lattice Constance of graphene, we will 
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have 

2^ ' ^^^^ 

where 6 is the argument of R. Therefore, we can conclude that the contact conductance is 
proportional to the square of the two contact areas between graphene and the two probes, 
and inversely proportional to the square of the probe interval, if the conditions given above 
are satisfied. 



4. Conclusions 



In this work we have systematically studied the contact conductance between graphene 
and two quantum wires which serve as the leads to connect electron reservoirs and graphene. 
The general conclusion we have obtained is that the contact conductance decays rapidly 
with the increase of the distance between the two leads. When each quantum wire couples 
to only one carbon atom in graphene, the contact conductance vanishes at the Dirac point 
if the two carbon atoms coupling respectively to the two leads belong to the same sublattice 
of graphene. And this conductance zero arises from the chirality of the Dirac electron in 
graphene. Also at the Dirac point, if two quantum wires couple to two carbon atoms of 
distinct kinds, a resonant path can be formed by adjusting the strength of the graphene-lead 
coupling. The inter-valley quantum interference causes the orientation dependence of the 
contact conductance. In the weak graphene-lead coupling limit, the carbon-carbon bond 
directions(i.e. Ti directions) are the optimal directions to form the maximal contact conduc- 
tance. In a more realistic situation, each quantum wire may couple effectively to multiple 
carbon atoms around it. To mimic such a situation, we introduce Gaussian-type graphene- 
lead coupling, by which we have worked out the two-dimensional conductance pattern formed 
by moving the second probe around the first one. We find that the conductance pattern does 
not vary sensitively with the incident electron energy. However, the symmetry of the conduc- 
tance pattern changes from C3 group to Cg group when the first probe shifts from a carbon 
atom to the center of a hexagon of the graphene lattice. Finally, we obtain an approximative 
expression about the contact conductance when the probe interval is sufficiently large and 
the incident electron energy is in the vicinity of the Dirac point. We find that in such a 
case the contact conductance is proportional to the square of the two contact areas between 
graphene and the probes, and inversely proportional to the square of the probe interval. 
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Appendix 

Now we derive the explicit expression about the self-energy terms in terms of the 
structure parameters. Substituting Eq. ffTOl) into Eq. ffTSj) we have 

W = ^, = ^v^iuF^^,{Ri^ - Rj^>)v*^,j^,. (35) 

In the above we have introduced four auxiliary functions F^^/ (R) which can be analytically 
treated by invoking the linear dispersion relation of graphene around the Dirac point. 

F^^{R)=F^^{R), (37) 
17 (T,^- ^ stk _VhM+ Ji{qR)q^dq 

and 



ks 



p 1^ e^^R stj _ V3^M_ J,{qR)q'dq 

where 7 = \/?>at/2 is the so-called Fermi velocity and 

M± = sin(^^±^). (40) 
3 

It should be noticed that in Eg. (1371) the incident electron energy E appears in the analytical 
result of the functions Fj^^{R) and Fgg{R) as a prefactor. It arises from the chirality of the 
Dirac electron in graphene. This result implies that these two functions are equal to zero 
when the incident electron energy is aligned with the Dirac point. Besides, our derivation 
indicates that the quantum interference between K and K' valleys in the band structure 
of graphene is responsible for the dependence of the four functions given above, hence the 
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contact conductance, on the relative orientation between the two probes. 
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FIG. 1: (Color online) (a) Schematic of the lattice structure of graphene and the reference frame. 
A hexagonal unit ceh represented by a dashed line contains two carbon atoms denoted by A and 
B. Three vectors directed from a A site to nearest neighbor B sites are given by ri = a[0, — 1/\/3], 
T2 = a[l/2, l/(2\/3)] and T3 = a[— 1/2, l/(2-v/3)] with a being the lattice constant. R is the 
distance between the centers of the two leads, (b) The first Brillouin zone of graphene. The two 
unequivalent vertices of the hexagon are called K and K' points, (c) An equivalent plot of the 
electronic transport structure shown in (a), but in (c) graphene is represented by two kinds of 
eigen-states of the vr band electron. 
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FIG. 2: (Color online) The linear conductance as a function of the incident electron energy for A-A 
coupling shown in (a) and (b), and A-B coupling shown in (c) and (d). In (a) and (c) the parameter 
£q takes several typical values which shift the energy band of the leads relative to the Dirac point 
of graphene. In (b) and (d) the parameter vq takes several typical values which correspond to 
different bandwidthes of the leads. 
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FIG. 3: (Color online) The linear conductance as a function of R, the interval between the two 
probes. The first probe is fixed at the origin. And the second probe departs from the origin along 
the y direction, (a) The conductance spectra in A-A(B) coupling for the incident electron energy 
E = 0.0, 0.1 and 0.5. (b) The case of Gaussian-type coupling for E' = 0, 0.1 and 0.5. 
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FIG. 4: (Color online) The linear conductance as a function of v, the graphene-lead coupling 
strength, (a) The case of A-B coupling, the incident electron energy is fixed at E = 0, but the 
probe interval varies from 11 = 2\ti\, |5ti| to |14ti|. (b) A-B coupling configuration with the probe 
interval fixed at -R = |5ti|, and the incident electron energy taking several different values: E=0, 0.1 
and 0.2. (c) Gaussian-type coupling configuration with the probe interval and the incident electron 
energy specified at E=0 and R = |5ti|. The decay factor in the Gaussian coupling function takes 
several different values: do = 0.4, 0.5, 0.8 and 1.2. 
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FIG. 5: (Color online) The two-dimensional conductance pattern formed by fixing the first probe 
at the origin labeled by a triangle, and shifting the second probe around, (a) The case of -E = 
and do = 0.2. (b) The case ol E = 0.5 and do = 0.2. 
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FIG. 6: (Color online) The two-dimensional conductance pattern formed by fixing the first probe 
at the origin labeled by a triangle, and shifting the second probe around, (a) The case of -E = 
and do = 1.0. (b) The case oi E = 0.5 and = 1.0. 
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